Численный анализ: от скаляра к вектору: вызов нелинейной системы

Переход от одного уравнения $f(x) = 0$ к системе с несколькими переменными — это ключ к решению сложных инженерных задач, от орбитальной механики до анализа структуры почвы. Мы больше не ищем простой ноль на линии, а ищем одновременное пересечение $n$ гиперповерхностей в $n$-мерном пространстве.

1. Математическая структура

Нелинейная система представляется как набор уравнений, где каждая компонентная функция зависит от вектора неизвестных $\mathbf{x} = (x_1, x_2, \dots, x_n)^t$:

$$f_1(x_1, x_2, \dots, x_n) = 0,$$ $$f_2(x_1, x_2, \dots, x_n) = 0,$$ $$\vdots$$ $$f_n(x_1, x_2, \dots, x_n) = 0,$$

Мы сводим это к векторной форме основная формула:

$$\mathbf{F}(\mathbf{x}) = \mathbf{0}$$

где $\mathbf{F} = (f_1, f_2, \dots, f_n)^t$. Индивидуальные функции $f_i$ обозначаются как координатные функции функции $\mathbf{F}$.

2. Аналитические основы и непрерывность

Для численного решения этих систем необходимо обеспечить хорошее поведение отображения. Определения 10.1–10.3 устанавливают, что пределы и непрерывность в $\mathbb{R}^n$ определяются компонентно.

Определение 10.3

Пусть $\mathbf{F}$ — функция из $D \subset \mathbb{R}^n$ в $\mathbb{R}^n$. Мы говорим, что $\lim_{\mathbf{x} \to \mathbf{x}_0} \mathbf{F}(\mathbf{x}) = \mathbf{L} = (L_1, L_2, \dots, L_n)^t$, если и только если:

$$\lim_{\mathbf{x} \to \mathbf{x}_0} f_i(\mathbf{x}) = L_i$$ для каждого $i = 1, \dots, n$.

Используя определение $\epsilon-\delta$: для любого $\epsilon > 0$ существует $\delta > 0$, такое что $\|\mathbf{F}(\mathbf{x}) - \mathbf{L}\| < \epsilon$ всякий раз, когда $0 < \|\mathbf{x} - \mathbf{x}_0\| < \delta$.

Ошибочная мысль: Независимость нормы

Критически важный нюанс: хотя можно использовать различные нормы ($\ell_1, \ell_2, \ell_\infty$), непрерывность не зависит от конкретного выбора. Существование предела не изменяется при любой векторной норме в $\mathbb{R}^n$.

3. Теоретический обзор

Теорема 1.6: Для функций из $\mathbb{R}$ в $\mathbb{R}$ непрерывность часто может быть показана путем доказательства дифференцируемости. В многомерном случае, если частные производные координатных функций существуют и ограничены, непрерывность гарантирована, что является необходимым условием для итеративных решателей.

Классический пример: Пример 1

Рассмотрим задачу о круглых пластинах в почве. Приведите систему из $3 \times 3$ нелинейных уравнений к стандартной форме $\mathbf{F}(\mathbf{x}) = \mathbf{0}$:

$3x_1 - \cos(x_2 x_3) - \frac{1}{2} = 0$
$x_1^2 - 81(x_2 + 0.1)^2 + \sin x_3 + 1.06 = 0$
$e^{-x_1 x_2} + 20x_3 + \frac{10\pi - 3}{3} = 0$

Здесь $\mathbf{x} = (x_1, x_2, x_3)^t$ и $\mathbf{F}(\mathbf{x}) = (f_1(\mathbf{x}), f_2(\mathbf{x}), f_3(\mathbf{x}))^t$.

ВОПРОС 1

Учитывая векторное отображение $\mathbf{F}(x_1, x_2, x_3) = (x_1 + 2x_3, x_1 \cos x_2, x_2^2 + x_3)^t$, почему $\mathbf{F}$ непрерывна в каждой точке $\mathbb{R}^3$?

Потому что это линейное отображение, и все линейные отображения непрерывны.

Потому что каждая координатная функция $f_i$ является суммой или произведением базовых непрерывных функций (полиномов и косинуса).

Потому что матрица Якоби имеет постоянный определитель.

Потому что она отображает $\mathbb{R}^3$ в замкнутое и ограниченное подмножество $\mathbb{R}^3$.

ВОПРОС 2

Какая из следующих функций $\mathbf{F}: \mathbb{R}^2 \to \mathbb{R}^2$ непрерывна всюду, кроме точки $(1, 0)$?

$\mathbf{F}(x_1, x_2) = (x_1 - 1, x_2)^2$

$\mathbf{F}(x_1, x_2) = (\sin(x_1 - 1), \cos x_2)$

$\mathbf{F}(x_1, x_2) = (1 / (x_1 - 1), x_2 / (x_1^2 + x_2^2))$

$\mathbf{F}(x_1, x_2) = (e^{x_1}, e^{x_2})$

ВОПРОС 3

Если мы меняем векторную норму с евклидовой нормы ($\ell_2$) на максимальную норму ($\ell_\infty$) при оценке непрерывности нелинейной системы, как это влияет на свойство сходимости последовательности векторов?

Последовательность теперь может расходиться, потому что максимальная норма менее чувствительна.

Сходимость сохраняется, потому что все векторные нормы в конечномерном пространстве эквивалентны.

Значение предела $L$ будет меняться в зависимости от масштабирования нормы.

Последовательность будет сходиться быстрее в норме $\ell_1$, чем в норме $\ell_2$.

ВОПРОС 4

Метод Ньютона обычно используется для нелинейных систем. Однако, если применить его к линейной системе $Ax = b$ (где $A$ невырожденная), какой результат получится?

Метод не работает, потому что якобиан не определён для линейных систем.

Метод сводится к методу Гаусса и требует $n$ итераций.

Метод достигает точного решения за ровно одну итерацию из любого начального приближения.

Метод колеблется и никогда не сходится, если $x_0$ не равно $b/n$.

ВОПРОС 5

При определении предела $\mathbf{F}(x)$ при $x$, стремящемся к $x_0$, что означает условие $0 < \|x - x_0\| < \delta$?

Мы оцениваем функцию точно в точке $x_0$.

Мы рассматриваем проколотую окрестность, то есть сама точка $x_0$ исключена из вычисления предела.

Расстояние от начала координат должно быть меньше $\delta$.

Система должна быть линейной для существования предела.